#include "PrimitiveCylinder.h"
#include "moremath.h"

Mat4d* PrimitiveCylinder::getMassTensor(){
	if(mMassTensor==NULL){
		mMassTensor = new Mat4d(vl_0);
		
		// WARNING - ellipsoid used for now : TODO
		(*mMassTensor)[0][0] = (mDim[0]*mDim[0])/5;
		(*mMassTensor)[1][1] = (mDim[1]*mDim[1])/5;
		(*mMassTensor)[2][2] = (mDim[2]*mDim[2])/5;
		(*mMassTensor)[3][3] = 1;
		(*mMassTensor)*=mMass;
	}
	return mMassTensor;
}
